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SUMMARY 


The paper considers dynamic plastic response of a thin cylindrical shell, 
immersed in a potential fluid initially at rest, subjected to internal 
pressure pulse of arbitrary shape and duration. The shell is assumed to re- 
spond as a rigld-perfectly plastic material while the fluid is taken as 
inviscid and incompressible. The fluid back pressure is incorporated into 
the equation of motion of the shell as an added mass term. Since arbitrary 
pulses can be reduced to equivalent rectangular pulses, the equation of mo- 
tion is solved only for a rectangular pulse. The influence of the fluid in 
reducing the final plastic deformation is demonstrated by a numerical example. 


FORMULATION OF THE PROBLEM 


Consider a rigid- ideally plastic, thin-walled, circular, cylindrical shell 
of Infinite length. The shell is surrounded by a pool of potential (inviscid 
and incompressible) fluid infinitely extended in all directions. The shell 
is subjected to an internal pressure pulse P(z,t), varying both along the axis 
and with time. P(z,t) is further assumed to be axisymmetric and S 3 mimetrlc in 
z with respect to z - 0 (fig. 1). 

This paper examines the influence of the fluid in reducing the plastic 
(residual) deformation of the shell. It is known that the pressure with 
which potential fluid resists the motion of a deforming solid can be con- 
sidered as an increase in the inertia of the solid. Therefore in order to 
solve the problem it is necessary to establish the so-called effective mass 
consisting of the actual mass of the shell and the added (virtual) mass 
reflecting the fluid resistance. Then the problem is reduced to the analysis 
of a shell deforming in vacuum. For the sake of continuity, we will adopt 
the notation introduced in reference 1. 


GOVERNING EQUATIONS 

The equation of motion of the shell is: 
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where M is the axial bending moment, N the circumferential (hoop) normal 
force, R, H and p the radius, the wall thickness and the mass density of the 
shell respectively, V the radial velocity of the points on the middle surface 
of the shell, and P(z,t) and Pf(z,t) are the externally applied pressure pulse 
and the back pressure of the fluid resisting the, motion of the shell 
respectively. 


We assume that the yield condition in the M,N space is defined by the 
limited interaction curve of fig. 2. The implications of this assumption are 
discussed in detail by Drucker (ref. 2) and Hodge (ref. 3). The yield values 
My and Ny are given by 

M = -y HRP N - RP (2) 

y 4 o y 0 

where 


P = 
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with Oy being the yield stress. 


(3) 


It is known (see, for example, ref. 1) that four different phases (modes) 
of plastic deformation may occur during the motion. We will consider herein 
only one of these phases which Occurs for all possible types of loading, 
though this restricts the magnitude of the loading to a certain limit. In 
the considered phase the deformation is characterized by a stationary plastic 
hinge circle at z = 0 and two moving hinge circles at z = + ?(t). 


It can be shown (see, for example, Eason and Shield (ref. 4)) that the 
plastic regimes (see fig. 2) are as follows; 


z = 0: M = -M 


z = C: M = M 


0 < z < C: -M < M < M , 

y y 


N * N 


N = N 


N = N 


Regime A 
Regime B 
Regime AB 


(4) 


Thus, from the normality of the strain-rate vector to the yield surface, 

2 “ 0 , V > 0 for 0 < z < 5 (5) 
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For this deformation mode the velocity, V(z,t), is therefore a linear function 
of z, i.e. , 


V(z,t) 



0 < z < Z 
z > Z 


( 6 ) 
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In the above equations and in the sequel because of symmetry it is enough to 
consider only half of the shell z ^0. 

DETERMINATION OF THE ADDED MASS 


Before attempting to solve equation (1), the backpressure Pf(z,t) should 
be determined as a function of V(z,t) and its derivatives. The equation 
governing the flow of the potential fluid is in polar coordinates 


.2 ^ r 9r . 2 ^ 


in 


r > R 


3r “ 9z 

where F(r,z,t) is the fluid velocity potential. 


(7) 


As the shell is impermeable, the velocities of the fluid and the shell 
at the points of contact must be identical, i-e., 

II = V at r = R (8) 

Furthermore, from the radiation principle, 

V -> 0 , - — >■ 0 , 0 as max(r, z) •> ” (9) 

or dz 

Once the fluid Velocity potential- is determined from the Laplace equa- 
tion (7), subject to the boundary conditions (eqs. (8) and (9)) the pressure 
exerted by the fluid on the shell can be computed from the Cauchy Integral, 

PfCz.t) = -p^ at r = R (10) 


where pf is the mass density of the fluid. 


The equations (7) and (10) imply the assxjmptlon that the perturbations 
about average values can be neglected. 

As a further approximation, we will assume that the functional relation 
between Pf(z,t) and t, is not sensitive to the time dependence of ?, and hence 
C may be treated as a constant for the determination of this relation. Then 
in view of equations (6), (7) and (8), we may write 


F(r,z,t) = V^(t) f(r,z) (11) 

and from equation (10), 

dV 

P^(z,t) = -P£f(R,z) (12) 

where -Pff (R,z) is the added (virtual) mass arising due to the resistance of 
the fluid being displaced by the shell. 

Substituting equations (11) and (6), (with ? being constant) into equa- 
tions (7) to (9), it follows 
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(13) 


9^f . 1 9f . 9^f 
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subject to 


li 

9r 


and 


■!) “ 

r = R, z 

< 4 

at 

II 

> t. 

9£ A 
— ->■ 0 , 

1^-0 

as 
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(14) 


(15) 


The details of the solution of equation (13) are omitted herein for the 
sake of brevity. A closed form integral solution is obtained after intro- 
ducing the Fourier cosine transform. The argument of this integral is 
rather complicated and the integration is performed in three stages using 
asymptotic formulae and Filon’s method, subject to the restriction, C < R 
which is subsequently seen to be not severe. In order to make this numerical 
solution amenable for substitution into equation ( 1 ), the result is subjected 
to a series of polynomial regression analyses. Finally we obtain 


where 

and 

with 


f(R, z) - -R <g 


0 VRy 1 \Rj 


gQ(x) = 3 q + 

2 3 

gj^(x) = Uq + a^x + + a^x 


Uq = .004994512 
= -.5420473 


“2 = 


,1058701 


= .1627119 


3q = . 02050149 
3 ^ = 1.664447 

82 = -1.105309 

83 = .4096600 


(16) 

(17) 

(18) 

(19) 


Note and 3^ are dimensionless constants that do not depend on the shell/ 
fluid parameters. 


PLASTIC DEFORMATION OF THE SHELL 


Equation (1) now becomes, in view of equations (2), 


= P(z,t) - P - 

9z^ 


PH + Pf«8olRjr\P«-Pf^l 

in 0 < z 


( 6 ), ( 12 ) and (16), 



< ? ( 20 ) 
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In equation (20), C — C(t). For arbitrary P(z,t) tbe above equation may only 
be solved by numerical methods. As a first step in simplifying the work, the 
approach introduced by Youngdahl (ref. 1) will be used to approximate a com- 
plex loading function by (i. e, , correlate it to) a simple rectangular pulse. 
Since the standard limit analysis of the shell is independent of any sur- 
rounding medium, the method given by Youngdahl (ref . 1) to determine the 
equivalent rectangular pulse does not need any modification in this case. 

Thus correlated, P(z,t) can be expressed as 


P(z,t) 

= P 

z < L and 0 < t < t 

e 

1 1 - g , 

F(z,t) 

= 0 

|z| > L or t > t^ 

' ' e e 


( 21 ) 


where Pg is the magnitude, tg the duration and 2 Le the length of the loaded 
area of the equivalent rectangular pulse (see ref . 1 for their derivation) . 


For plastic deformation to occur, Pe must be greater than the limit load. 
This condition is expressed by (see ref. 1) 


P 

e 



( 22 ) 


For the deformation to take place in the assumed phase, the following 
boundary conditions must be satisfied 


M = -M , ^ = 0 at z = 0 

y oz 

3M 

M = My , If = 0 at z = ? 


(23) 


Further, the condition that the bending moment does not exceed My at the 
hinge circle at z = 5 implies 


^ < 0 at z = C ( 24 ) 

3z 


In addition the condition that the bending moment cannot be less than -My at 
the hinge circle at z = 0 Implies 



> 0 


at z = 0 


(25) 


For the interval 0 < t < t^, a trial solution as in the case of a shell 
deforming in vacuum is assumed. This is taken in the form 


md 


c(t) - 2 ^ 


V^(t) 



0 i t < t 

e 


(26) 


fhere z^ and are constants. Substituting equation (26) into equation (20) 
ind integrating twice subject to the boundary conditions (23) yields in the 
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end two equations for zj[ and These two can be reduced to; 


and 


where 


a(zJP z? + {6 - 4a(z,)}P L z, - {6 - 3a(z ') }(P L^ + P BH) 
1 o 1 1 eel lee o 


= 0 


^1 = 


6pH(P L z, -P L - P RH) 
eel ee o 


2 

zi PH-PfRgi 


a(Zj^) 


fh] 

pH P£ ^^1 1 R ) 


PH + Pj Rg 


0 (rJ 


The inequalities (24) and (25) can be written in the form 

Vi + T i (*o (t) ^ (^)) 4 ’ ® 

and 4P L z. - P z^ - 3P - 3P„RH < 0 

eel el e e ° 


(27) 


(28) 


(29) 


(30) 

(31) 


From equations (17) , (18) and (19) it is easily verified that inequality 
(30) is always satisfied. If inequality (31) is not satisfied motion cannot 
start in the phase assumed. When the inequality becomes an equality, Pg takes 
the bounding value. To determine the bounding value of Pg, the non-linear 
equations (27) and (31) should be solved simultaneously. This is done 
numerically. Figure 3 shows the range of values of P^ that satisfies in- 
equalities (22) and (31) and hence gives rise to deformation in the assumed 
mode. For a Pe belonging to this range, the non-linear equation (27) may be 
solved numerically. Also it may be verified that the restriction z^ < R is 
always satisfied if Lg < R. Thus, the solution discussed in this paper is 
valid for Lg < R and Pg satisfying inequalities (22) and (31) . 


For t > tg, there is no internal pressure. Letting P(z,t) = 0 in equa- 
tion (20) , integrating it twice with respect to z and substituting the result 
into the boundary conditions (23) , we arrive at the following equations 


dVo -P (?2 + 3RH) 

_ = ___ = 


(32) 


J_ 

dt 



-6PqRH 

5 l^pH- p^Rg^ 


4>a) 


(33) 


Equations (32) and (33) constitute a system of non-linear, first order differ- 
ential equations for and Z, the initial conditions being given by 
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(34) 


C(tp = 

K,t 

V (t ) = 
o'^ pH 


The above differential equations are valid only for t^ s t < t^, where is 
defined by 

V (tp) = 0 (36) 

o f 


We will denote 

= ?(tj) 

From equations (32) and (33), we can express as, 

V (ti . 3M+ii_\ 

° ^ (fl) (,1) 

From (36) and (37), we see that can be obtained from the equation, 

3RH + 


PH-P^ 1-^ 


pH + p^ Rgg I ^ 


Equation (39) can be solved numerically to obtain ^f. It is noted that 
Cf depends only on the shell parameters H and R and the density ratio pf/p. 
Since Cf is the quantity that is known and not tf, the equations (32) and 
(33) are now reformulated with ^ being the independent variable and t(^) and 


V (5) being the dependent variables. Thus, 
o 


d? dC 


dC I ? 


o \ dt 


Zi < 5 ^ 


The new initial conditions are 




VJzp - 


K,t 
1 e 


Equations (40) and (41) are solved numerically. Finally, the maximum 
plastic deformation U^(t) can be obtained as. 


' -w 


0 < t < t 
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ft 


and 


U (t) = U (t ) + 
o o' e 


VQ(t)dt 


0 < t < t. 


or 


Ki 

= -U ^ 


V^(S) ^ d? 

di 


Zi < ? < Sf 


( 45 ) 


The integral in equation (45) can be nximerically evaluated after the numeri- 
cal solutions V^(?) and t(c) have been obtained. 

For the special case in which zi coincides with Cf , the solution for 
t > tg discussed above is not valid. For this case, equations (32) and (33) 
with (35) yield. 


= ' "i ■ 'f 

Vo(t) - 


(46) 


(47) 


where <j>(C) is defined in equation (32). Note < 0. From equation (47) 

we see 


K, 


^f = "e 1 


PH<|) (z^) 


(48) 


From equations (47) and (44) we can show that 

^(z.) 2 ^l^e 

Uo(t)=-^(t-y +— (2t-y 


(49) 


Finally we have, for the maximum plastic deformation in this special case 

^ - K. 


2pH " pH$(z,) 


(50) 


NUMERICAL EXAMPLE 


For a shell with H/R = 1/36, Lg/R =1/4 surrounded by a fluid of 
pf/p = 1/10, the complete numerical solution is determined for the admissible 
range of loads Pg. As is seen from figure 3, the range of Pe/Po that will 
give rise to motion in the assumed phase is between; 1.33 and 2.97. The 
same range for a shell in vacuum is 1.33 to 2.19. Figure 4 shows the final 
maximum plastic deformation, Uo(tf) , (non-dimensionalized as pHUo/Fo*^e ) ^ 

function of Pg/Po* For the sake of comparison the corresponding curve for a 
shell deforming in vacuum is also shown in the same figure. 

In the numerical methods used, non-linear algebraic equations such as 
(27) and (39) were solved by Newton’s iteration method and the system of 
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differential equations (40) , (41) by a method using automatic step change 
(ref. 5). 
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Figure 3.- Range of pulse intensity 
initiating motion in assumed phase 
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Figure 4.- Maximum plastic deformation 
as a function of pulse intensity. 
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